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Re-localization due to finite response times in a nonlinear 
Anderson chain 
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Abstract. We study a disordered nonlinear Schrodinger equation with an additional relaxation process 
having a finite response time r. Without the relaxation term, r = 0, this model has been widely studied 
in the past and numerical simulations showed subdiffusive spreading of initially localized excitations. 
However, recently Caetano et al. (EPJ. B 80, 2011) found that by introducing a response time r > 0, 
spreading is suppressed and any initially localized excitation will remain localized. Here, we explain the 
lack of subdiffusive spreading for r > by numerically analyzing the energy evolution. We find that in the 
presence of a relaxation process the energy drifts towards the band edge, which enforces the population of 
fewer and fewer localized modes and hence leads to re-localization. The explanation presented here relies 
on former findings by the authors et al. (PRE 80, 2009) on the energy dependence of thermalized states. 
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Localization in linear, disordered systems is a well es- 
tablished phenomenon also known as Anderson localiza- 
tion [3]. One interesting aspect is the interplay between 
localization due to disorder, and nonlinearity, which by in- 
tuition should lead to thermalization in high dimensional 
systems. This is mostly formulated as a spreading problem 
of initially localized states and has first been addressed 
mainly numerically [4l, but lately also experimental stud- 
ies ^SjJ)^ as well as mathematical treatments [3^ were estab- 
lished. Numerical results showed a sub-diffusive spreading 
of an initially localized wave packet over a large timescale 
up to i ^ 10^ [Slin]- Lately, the spreading behavior was 
also related to properties of high-dimensional chaos fTO^ . 

There are a number of models where disorder and non- 
linearity can be studied, for example Hamiltonian oscilla- 
tor chains [11], the Klein-Gordon model [8| or coupled map 
lattices [12]. However, here we use the Discrete Anderson 
Nonlinear Schrodinger Equation (DANSE) and we start 
with a short summary of its properties found in previous 
studies [2]- The equations of motion for this model are: 



ilpn = V„1p„ + 1pn-i -t- -0„+i -I- fi\lpn\^ll}n 

and can be derived from the Hamiltonian: 



(1) 



H 



? + M*n+l + + ^IV'nl"- (2) 



■(/;„ is the complex valued field amplitude, 1^ is a random 
potential choosen iid. from {—W/2 . . . W/2) with typically 
W — 4.0 and /3 describes the nonlinear strength with a 
typical value of /? — 1. We emphasize that this model 
exhibits two conserved quantities: norm A/" = ^ iV'nP = 1 
and energy E = H . While the norm can be scaled to 



= 1 by a rescaling of /?, the energy is an important 
quantity governing the behavior of the system. This en- 
ergy dependence has been investigated in [2] and we shortly 
summarize the results, as they will help understanding the 
modified model studied below. In the linear case, /3 = 0, 
the Hamiltonian ([2]) can be diagonalized: H = e^c^Cfc 
whith energy eigenvalues efc and eigenfunctions ipk where 
ipn = X]fc Ck^k,n- Due to the disorder 14,, Cfc are randomly 
distributed in the interval [—W/2 — 2, W/2 + 2] and the 
are exponentially localized in space |3] (for a review see 
e.g. [13] )• Adding nonlinearity /? 7^ one still can write the 
system in terms of the linear eigenfunctions, but now an 
additional 4-mode coupling term will appear in the diag- 
onalized Hamiltonian allowing for energy / norm exchange 
between linear modes, hence spreading. 

In ^ it was shown that the energy is a crucial parame- 
ter determining the efficiency of spreading. If one initially 
excites an eigenmode with energy at the band edge, no 
spreading can appear because of energy and norm con- 
servation. There simply does not exist a state with such a 
norm and energy value other than the localized eigenmode 
at the band edge. If, on the other hand, time evolution is 
started from an excitation with an energy close to the 
band center i? « 0, the conservation laws do not prevent 
the system from reaching extended states and typically 
spreading is observed numerically. This can be formulated 
more rigorously by considering thermalization in such a 
system in terms of a Gibbs distribution [2]. It should 
be noted that a thermalized state is usually also maxi- 
mally spread as a maximum Gibbs-entropy corresponds 
to a maximum number of excited lattices sites or eigend- 
modes, respectively. Denoting Ck the amplitude of the fc-th 
eigenmode such that ipn — <^kfk,n we find for the ther- 



2 



Mario Mulansky, Arkady Pikovsky: Re-localization due to finite response times in a nonlinear Anderson chain 




time time 



Fig. 1. Energy of the state as function of time for /3 = 1.0 with r — 0.1 (left) and r = 0.001 (right). Shown is the time 
evolution of an initially thermalized state (maximal entropy, Gibbs distribution, see text) in a lattice with A'^ = 32 sites for a 
given temperature T = —100 . . . 100. 



malized state: 

\ck\^ = e-'''/^/Z, with Z = J2e-"'/^, (3) 

k 

where T £ (—00,00) is the dimensionless effective temper- 
ature. Note that for a specific disorder realization, fixing 
T is equivalent to fixing the energy E, so that tempera- 
ture values ±00 correspond to energies at the band cen- 
ter while T = ±0 belong to energies at the band edges. 
From ^ it is seen that in the limit T — > ±00 the ther- 
malized state is an uniform excitation of all eigenmodes, 
hence extended, while for T ±0 only the eigenmode 
with largest/smallest eigenvalue gets excited which results 
in a localized final state. So changing the energy from the 
band center towards the band edge makes the thermalized 
state becoming more and more localized. This has been 
found in [2J and will now be used to explain the results in 
a modified model. 

Recently, Caetano et al. introduced a modification of 
the DANSE model that consists of an additional relax- 
ation process introduced on top of the nonlinearity Ij: 

ii>„ = Vn^Pn + i>n-l + + P Mn 

Tfn = IV'nP - fn- 

It was found numerically by Caetano et al. that in this 
system, any non-zero relaxation time r > surpresses 
spreading. However, up to now no explanation for this 
absence of spreading has been given. We now are going to 
present such an explanation. 

We first note that by introducing the relaxation pro- 
cess we destroy the energy conservation valid for the orig- 
inal DANSE model (t = , /„(t) = However, 
the norm A/" = J2 IV'nP is still conserved here and we again 
fix A/" = 1. The energy becomes time dependent and cal- 
culates as: 

m = E + m^nrn+l) + f /^Wl^nl'- (5) 

n 



Numerical simulations presented below show that for 
any r > the energy increases with time if /3 > 0, while 
it decreases for /3 < oQ But as M is still conserved in this 
system, the energy can not increase (decrease) unbounded. 
Norm conservation still enforces a strict bound on the en- 
ergy E, namely -2 - W/2 + P/2 < E <2 + W/2 + P/2. 

In a first numerical study we fixed a lattice with = 32 
sites and used a thermalized state as given by ([3]) with a 
given temperature T as an initial condition for il)n and 
/n(0) = 0. Starting from such an initial condition for a 
given temperature and fixing (3 — -1-1, we integrated the 
equations of motion ^ using the DopriS-Runge-Kutta 
method of order 5 with error control (tolerance 10~*) and 
dense output [M] . This is not the most efficient numerical 
method for a stiff system like our modfied DANSE, but 
for the results presented here the computational effort did 
not require us to apply more sophisticated routines. Ex- 
emplary results for t = 0.1 and t = 0.001 are shown in 
Fig. m For all tested values of r = 0.001 ... 100 we found 
an energy growth during the time evolution up to a satu- 
ration value Eg ~ 3.0, which corresponds to the maximal 
energy eigenvalue, hence the upper band edge. Obviously, 
the time-scale at which this energy drift takes place is gov- 
erned by T. As seen in Fig. [TJ changing r by two orders of 
magnitude also changes the time-scale of the energy drift 
by roughly two orders of magnitude. Numerical results 
for P — —\ are very similar but now the energy decreases 
down to i?s ~ —3 for all values of r and all tested initial 
conditions. 

From the consideration above we can now conclude 
that even if thermalization occurs, the final state must 
be localized as it will always have an energy close to the 
band edge. So the numerical observation of the drift of the 
energy together with the previously found energy depen- 
dence of thermalized states [2] explains the lack of spread- 
ing in systems with finite response time. We note that the 



^ Note that we used /n(0) = in all simulations, for other 
initial values /n(0) the energy behavior might also depend on 
the sign of the term \ijjn\^ — fn- 
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Fig. 2. Spreading ol initial eigenstates with lowest (/J = 1) and highest (/3 = —1) energy value for r = 0.01 (left panel). In the 
right panel we initialized the system with eigenstates from the band center (E ~ 0) and started time evolution with P — ±1 
using T — 0.0001. Shown are participation number P ~ i/Yl IV"".!* s^nd energy as functions of time. 



localization of high (low) energy modes can also be under- 
stood in terms of a self-trapping mechanism. Indeed, if the 
energy is shifted out of the spectrum of energy eigenvalues 
of the linear problem (/3 = 0), then by adding nonlinearity 
a localized mode is created and spreading is prevented [T^ . 

In a second numerical experiment we revisit the spread- 
ing problem for system ^ to further visualize the inter- 
play between spreading and energy drift. We took a lat- 
tice of = 256 sites and choose either the eigenmode with 
lowest /highest energy or an eigenmode from the band cen- 
ter as initial condition for ipn and again /n(0) = 0. Then 
we ran the time evolution with /3 = ±1.0 and measured the 
spreading as well as the energy E. As measure of spread- 
ing we chose the participation number P = l/X)!''/'"!^ 
that is a widely used quantity in such problems as it 
roughly counts the number of excited sites. Large num- 
bers of P correspond to spreading while small P indicate 
localization. The results are shown in Fig. [21 As is seen 
there, starting from the highly localized eigenstates at the 
band edge the relaxation process increases (/3 — +1) or 
decreases (/3 = —1) the energy which initially leads to 
spreading. When the energy crosses zero the wavefunction 
is maximally delocalized, but a further increase/decrease 
leads back to localization until the energy reaches the 
other edge of the energy band and the process saturates 
in a localized state. Hence, the energy drift towards the 
band edges induces re-localization of the wave function. 
The same happens when starting from states with initial 
energies at the band center i? « as also shown in Fig. [51 
There, some initial spreading takes place due to nonlinear 
interactions, but finally also re-localization appears due to 
the energy drift. 

To summarize, we found an explanation for the sur- 
pression of subdiffusive behavior in a modified DANSE 
model with an additional relaxation process with response 
time r > 0. Our reason is based on the numerical observa- 
tion of an energy drift induced by the relaxation process. 
This energy drift continues up (down) to the energy bound 
as given by norm conservation in the system. Using previ- 



ous results about thermalization properties in the DANSE 
model stating that states close to the band edges must be 
localized |^ , we found the explanation for the effect of re- 
localization observed here and in a previous study [1] . We 
note that our explanation is based on the numerical ob- 
servation of an energy drift. A rigorous analytic treatment 
of the dynamic equations should be done to obtain math- 
ematical results showing the drift of the energy, but this 
is left for future work. To our opinion, the importance of 
the energy on thermalization as described in [2] is crucial 
for understanding such modified DANSE models. 
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